Ga substitution as an effective variation of Mn-Tb coupling in multiferroic TbMnOa 
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Ga for Mn substitution in multiferroic TbMnOs has been performed in order to study the in- 
fluence of Mn-magnetic ordering on the Tb-magnetic sublattice. Complete characterization of 
TbMni_a;Gaa;03 {x = 0, 0.04, 0.1) samples, including magnetization, impedance spectroscopy, and 
x-ray resonant scattering and neutron diffraction on powder and single crystals has been carried 
out. We found that keeping the same crystal structure for all compositions, Ga for Mn substitution 
leads to the linear decrease of Tj^" and r^", reflecting the reduction of the exchange interactions 
strength Jmii-Mh and the change of the Mn-O-Mn bond angles. At the same time, a strong suppres- 
sion of both the induced and the separate Tb-magnetic ordering has been observed. This behavior 
unambiguously prove that the exchange fields Jivin-Tb have a strong influence on the Tb-magnetic 
ordering in the full temperature range below Tj^" and actually stabilize the Tb-magnetic ground 
state. 

PACS numbers: 75.47.Lx, 75.25. +z, 75.30.Et, 75.30.Kz, 75.80.-|-q 



I. INTRODUCTION 

Rare earth manganites i?Mn03 (i? ~ Gd, Tb, Dy, 
Ho) with orthorhombically distorted perovskite structure 
have been attracting a lot of attention due to a variety 
of unusual physical properties including a potentially in- 
teresting cross coupling of magnetism and ferroelectricity 
(FFi)JuL^ These so called multiferroics are of particu- 
lar interest for understanding the fundamental physical 
links between spin, charge and lattice degrees of freedom 
that give rise to magnetoelectric coupling, as well as be- 
cause of the promising possibility of using these coupled 
order parameters in novel device applications by control- 
ling the material's polarization state with either electric 
or magnetic fieldi^i^ii 

While it has been clearly shown, that the ferroelectric- 
ity in these compounds is magnetically driven, the role 
of magnetic 3d and 4/ ions was found to be rather differ- 
ent. As commonly accepted, the origin of multiferroicity 
in i?Mn03 is a result of a cycloidal Mn-magnctic order- 
ing with inverse Dzyaloshinsky-Moriya interaction being 
the driving force of polar lattice distortions i^'^'^^i^^'^^i^^ 
Spontaneous electric polarization can be therefore ex- 
pressed via the niy and rriz components of the Mn-spins 
of the cycloid and the magnetic propagation vector, r, 
as Fs c)c mymz{ex x r), where ex is the unit vector 
along the axis of rotation of Mn-spinsjiS Magnetic rare 
earths in turn, may determine the polarization direc- 
tion via their interaction with Mn-spins. In compounds 
with rare earths showing strong magnetic anisotropy, the 
polarization is oriented along c-axis, P||c for R = Tb 
and Dy, while in those with nonmagnetic (Eu,Y) or less 
anisotropic Gd, P|| ai^ii^ In some cases, they even con- 
tribute to the magnitude of the polarization as has been 
shown for DyMnOsj^ii^ A necessary condition for this 
is the ordering of rare earths with the same propagation 



vector as Mn, behavior confirmed now for i?Mn03 with R 
= Gd, Tb, Dy and Ho above their own ordering temper- 
atures, < 10 K,^J^iiI Here, Mn -magnetic sublattice 
polarizes i?-spins and the induced i?-moment is defined 
by the strength of exchange interaction between Mn- and 
i?-ions, Jmii-r- 

As proposed in Rei^ the ordering of the i?-spins be- 
low (often referred as individual) is strongly depen- 
dent on the relative strengths of exchange interactions be- 
tween Mn- and i?-ions, Jmh-r, and between rare earths 
themselves, Jr_r. will not significantly affect the R- 
ordering as found for R = Dy that orders on its own 



^ r'^") below 6 K,^^ On the other hand. 



strong Jmii-r will force the i?-ordering to the same peri- 
odicity as Mn down to the lowest temperatures as in case 
of ii;= Ho (rH° ^r"")^^ 

The most interesting case, however, is an intermedi- 
ate coupling regime as observed in TbMnOs. Here, the 
Tb- and Mn- orderings remain coupled down to the low- 
est temperatures through the harmonic coupling of their 
wave vectors, ir'^^ — t'^''" = 1. This is a result of ad- 
justment of Ising-like Tb spins to periodic Mn ordering 
minimizing system's energy and leading to rather com- 
plex Tb-spin structure at low temperatures This inter- 
mediate coupling regime should be very sensitive to any 
variation of Jmh-r and, correspondingly, strong effects 
on the Tb-magnetic ordering are expectedi^ 

In this paper we report on combined neutron diffrac- 
tion, x-ray resonant magnetic scattering, and single crys- 
tal magnetization and dielectric measurements on sub- 
stituted TbMni_3;Ga^03 {x = 0, 0.04, 0.1) compounds. 
Substitution of Mn^+ ion by the closest isovalent but 
nonmagnetic Ga'^+ was chosen to provide minimal distur- 
bance of the lattice and effectively vary the total Jmh-Mh 
and; consequently, Jmh-r, while keeping Jr_r constant. 
We show that; while keeping the same crystal structure 
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for all compositions, Ga for Mn substitution leads to the 
linear decrease of Tj^" and r'^" reflecting the intended 
decrease of Jmh-Mii and the increase of the Mn-O-Mn 
bond angles. At the same time we observe a strong sup- 
pression (for X = 0.04) and disappearance (for x > 0.1) 
of both the induced and the individual Tb-magnetic or- 
dering. This behavior confirms that the exchange fields 
^Mn-Tb from the Mn have a strong influence on the Tb- 
magnetic ordering in the full temperature range below 
the ferroelectric transition and actually stabilize the Tb- 
magnetic ground state. 



II. EXPERIMENT 

Polycrystalline TbMni-^Ga^rOa {x = 0, 0.04, 0.1) 
samples were prepared from a mixture of Tb407 (4N), 
Mn203 (3N) and Ga203 (4N) using standard solid state 
reaction. The lowest content of Ga {x i2 4%) was cho- 
sen such as to obtain a homogeneous distribution of Ga 
within the sample. The highest content {x ±2 10%) was 
determined by our ability to obtain a single phase sam- 
ple at this composition without significant changes in its 
crystal structure and magnetic properties. The single 
crystals were grown in argon atmosphere by the floating 
zone technique in the IR-heated image furnace (NEC) 
equipped with two halogen lamps (500W) and double el- 
lipsoidal mirrors. The growth and the feed rates were 
maintained to obtain a stable molten zone. 

Crystal structure and magnetic ordering of 
TbMni_^Ga3;03 were investigated between 2-300 K by 
neutron diffraction (ND) on powder and single crystal 
samples using the high resolution E9 diffractometer 
(A =1.797 A and 2.816 A) and the El thermal triple axis 
spectrometer operated in two-axis mode (A =2.44 A), 
respectively, at the Helmholtz-Zcntrum-Bcrlin (HZB). 
The data were analyzed with the FuUProf refinement 
package.— Single crystal x-ray magnetic resonant scat- 
tering (XRMS) on TbMni_2;Gax03 were conducted 
at the 7 T multipolc wiggler beamline MAGS, also 
HZB^^ 

Physical characterization of the samples was performed 
on both single crystal and polycrystalline samples. The 
single crystals were aligned along their principal direc- 
tion by x-ray Laue diffraction and cut in orthorhombic 
shapes with each face normal (within 2°) to a principal 
crystallographic direction. 

Temperature dependence of low and high field magne- 
tization along three main crystallographic directions were 
measured on TbMni_a;Gaa;03 {x = 0, 0.04, 0.1) in the 
VSM mode of a Physical Properties Measurement System 
(Quantum Design PPMS). The dielectric measurements 
were performed in a 14.5 T Oxford Instruments cryomag- 
net equipped with a ^He Heliox insert. Opposite faces of 
the plate-like samples were sputtered with gold and lo- 
cated between two gold covered copper plates. The top 
plate and the sample were pressed to the fixed and ther- 
malized bottom plate by a soft spring. The capacitance 



of such a plate capacitor with a sample as a dielectric 
was measured by an Andeen-Hagerling 2500 A capaci- 
tance bridge working at a frequency of 1 and 20 kHz. 

III. RESULTS AND DISCUSSION 

A. Crystal structure 

Starting with the structural analyses, powder ND 
measurements at room temperature show that all 
TbMni_:cGa3;03 samples are single phase and crystal- 
lize with the orthorhombically distorted perovskite struc- 
ture (space group Vhnm). Ga shares Ah atomic position 
(Wyckoff notation) with Mn and its refined content is 
in good agreement with the nominal one (Table We 
would like to emphasize that despite to close scattering 
lengths of Tb (7.38 fm) and Ga (7.29 fm), sharp differ- 
ence between those of Mn (-3.73 fm) and Ga allows to 
determine their occupancies confirming that all Ga is ac- 
counted in the Mn-site. With increasing Ga-content up 
to 10% (Fig. [T^-b), the lattice linearly contracts along 
fo-direction, Ah/b ~ 0.7%, and expands along the c-axis, 
Ac/c ^ 0.3% while leaving the a-parameter almost con- 
stant, Aa/a < 0.1%. The resultant decrease in volume 
is about 0.4%, approaching the volume of heavier R - 
DyMnOs. Anisotropic change of the lattice, however, 
leads to an increase of Mn-O-Mn bond angle (by ^ 0.4°) 
as if moving the system to the lighter R - GdMn03. 
Disturbance of the original TbMnOs lattice is also re- 
flected in changes of Mn-0 and Tb-0 interatomic dis- 
tances, both having maximal Ad/d 1.1% as presented 
in Table ni 

B. Magnetic and dielectric properties 

Turning to the magnetic properties, in Fig. [2] we 
present the temperature dependence of the magnetiza- 
tion of TbMni_2;Ga2;03 single crystals measured along 
the a- and 6-crystallographic directions with H = 5 kOe 
and along c with 1 kOc. We start with magnetization 
curves measured with H\\c as they show anomalies cor- 
responding to the onset of both Mn- and Tb-magnetic 
orders. As previously reported for the a; = compound. 
Ale presents three sharp anomalies at Tj^" = 42 K, 
Tfe = 27 K and T^^= 7 K attributed to the Mn Neel 
transition temperature, the transition temperature to 
the Mn-cycloidal state associated with development of 
ferroelectricity and the Tb Nccl transition temperature, 
respectively!^ The 4% Ga-containing compound preserves 
all three transitions. The first two are shifted down to 
39 K and 25 K whereas the last anomaly is enormously 
reduced and has a maximum around 4.5 K. The 10%-Ga 
sample shows only one sharp transition at T^" = 36 K 
and broad anomaly at around 15 K which might reflect 
the remanent ferroelectric transition TpE- Any signature 
of T^^ is clearly absent at this composition. 
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Temperature 




300 K 






12 K 




X 





0.04 


0.10 





0.04 


0.10 


a (A) 


5.30159(4) 


5.30049(7) 


5.29981(7) 


5.31606(4) 


5.31406(6) 


5.31268(7) 


b(A) 


5.85557(4) 


5.83318(7) 


5.81546(7) 


5.82304(5) 


5.80384(8) 


5.78569(7) 


c(A) 


7.40003(6) 


7.41209(9) 


7.42329(9) 


7.38443(7) 


7.39939(8) 


7.41051(11) 


x(Tb) 


-0.0150(3) 


-0.0160(4) 


-0.0162(4) 


-0.0162(4) 


-0.0159(4) 


-0.017(4) 


y(Tb) 


0.0814(3) 


0.0801(3) 


0.0796(4) 


0.0789(3) 


0.0790(3) 


0.0778(4) 


x(Oi) 


0.1055(4) 


0.1064(5) 


0.1050(5) 


0.1064(4) 


0.1058(4) 


0.1059(5) 


y(Oi) 


0.4658(4) 


0.4673(4) 


0.4667(5) 


0.4666(4) 


0.4686(4) 


0.4678(4) 


X(02) 


0.7035(3) 


0.7038(4) 


0.7029(4) 


0.7037(3) 


0.7038(3) 


0.7046(4) 




U.OZ i\jyo) 


U.OZD ( yo ) 










Z(02) 


0.0515(2) 


0.0514(2) 


0.0515(2) 


0.0515(2) 


0.0519(2) 


0.0519(2) 


Occ(Ga) 





0.044(2) 


0.104(6) 





0.044 


0.104 










0.2754(2) 


0.2651(3) 


0.2463(5) 


mf- 








3.95(4) 


3.30(4) 


3.19(8) 


mf" 








2.46(10) 


1.86(15) 


0(1) 


t"""" (b*) / 2 K 








3/7 






m^^ / 2 K 








5.79(6) 






m^V 2 K 








3.11(8) 






mJV 2 K 

















4.67 


3.83 


3.79 


4.14 


3.26 


5.47 


r>Mn 
^mag 








6.99 


6.47 


7.0 


pTb 
-*^mag 








13.7 








1.97 


3.17 


2.62 


1.69 


3.49 


3.33 



TABLE I: Crystal and magnetic structure parameters for TbMni_a;Gaa;03 (a- = 0, 0.04, 0.1) obtained from the Rietveld 
refinements of powder ND at room and low temperatures using 1.797 A and 2.816 A. The atom positions are given for Phnm 
settings where Mn occupies 4& (1/2 0), Tb and Oi 4c {x y 1/4), and O2 8d {x y z) Wyckoff positions. The low temperature 
data for all samples are given for 12 K except of Tb-magnetic structure parameters given for 2 K. It is caused by difficulties in 
fitting the Mn-magnetic refiections being overlapped with Tb-magnetic peaks below 8 K. 



X 





0.04 


0.10 


Mn-Oi (A) 


1.9431(7) 


1.9463(8) 


1.9471(8) 


Mn-Oa (A) 


1.9069(16) 


1.9060(20) 


1.9131(20) 


Mn-02 (A) 


2.2331(17) 


2.2235(18) 


2.2084(18) 


Mn02Mn (°) 


144.99(7) 


145.15(8) 


145.22(8) 


Mnd^4n (°) 


144.39(3) 


144.38(3) 


144.77(3) 


Tb-Oi (A) 


3.661(3) 


3.633(3) 


3.622(4) 


Tb-Oi (A) 


2.340(3) 


2.350(3) 


2.341(4) 


Tb-Oi (A) 


3.203(3) 


3.198(3) 


3.189(3) 


Tb-Oi (A) 


2.274(3) 


2.269(3) 


2.276(3) 


Tb-02 (A) 


2.542(2) 


2.538(2) 


2.535(3) 


Tb-02 (A) 


2.570(2) 


2.578(2) 


2.582(2) 


Tb-02 (A) 


2.317(2) 


2.311(2) 


2.311(3) 


Tb-02 (A) 


3.666(2) 


3.655(2) 


3.648(3) 



TABLE II: Selected distances and angles from Rietveld refine- 
ments of powder ND data for TbMni-^Ga^Oa {x = 0, 0.04, 
0.1) at RT. 



Magnetization curves along a- and 6-axis for all compo- 
sitions look very similar showing mainly a paramagnetic 
Curie- Weiss behavior of Tb-moments down to T^'^ as 
shown in the inserts in Fig. [21 The anomaly correspond- 
ing to onset of Tb-magnetic ordering at T^^ shifts to 
lower temperatures and becomes broader with increasing 
Ga content. 

Having observed a rather similar magnetic behavior 
of the Mn-sublattice in the pure and the substituted 
TbMnOa compounds, we have performed dielectric mea- 
surements for the X = 4% and 10% samples. Impedance 
spectroscopy measurements from 1 kHz up to 20 kHz re- 
veal no relaxor behavior for all compounds. Similar to 
TbMnOa, the 4% and 10% Ga samples present a strong 
anisotropy and anomaly in the dielectric constant along 
c- and a-directions at 40 K. In TbMnOs the ferroelectric 
transition is characterized by a lamda-shaped anomaly in 
the c-direction at 27 K (Fig.[3l).-^ The dielectric measure- 
ments along the c-direction for 4% and 10% samples show 
similar anomalies shifted to lower temperatures. They 
appear at 24 K and 15 K, respectively, in perfect agree- 
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FIG. 1: (Color online) Magneto-structural parameters of 
TbMni_a;Gaa;03 compounds as determined from ND: lattice 
parameters and volume (a)-(b) at RT, magnetic propagation 
vector at 2 K, Mn Neel temperature and ferroelectric transi- 
tion temperature (c), FWHM of Tb magnetic reflections and 
correlation length of Tb clusters at 2 K (d) as a function of re- 
fined Ga-content. Coherence length has been estimated using 
the Scherrer formula. Cross symbols in panel (d) correspond 
to instrument resolution at the position of the strongest Tb- 
magnetic reflection 
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FIG. 2: a) Temperature dependence of magnetic d.c. sus- 
ceptibility for TbMni_iGaj;03 {x = 0, 0.04 and 0.1) single 
crystals measured at 1 kOe along c-direction. The insert rep- 
resents the d.c. susceptibility of TbMni-ajGa^Os single crys- 
tals measured along the a- and &-crystallographic direction 
with H — 5 kOe. Data have been collected on heating after 
previous zero field cooling down to T=2 K. 



ment with TpE determined from magnetization measure- 
ments. The height of the lamda-shaped anomaly of the 
dielectric constant {6e) is reduced with increasing Ga con- 
tent from 5e of about 4 for x = down to 1.8 and 1 for 
X =4% and 10%, respectively. All these features suggest 
that at the respective TpE the compounds undergo a fer- 
roelectric transition that becomes less pronounced with 
increasing Ga-doping. 

The results from the measurements of magnetization 



FIG. 3: Temperature evolution of relative dielectric constant 
(e^) for TbMni-^cGa^Os with a; = 0, 0.04, 0.1, respectively, 
measured on cooling with E — 20 kHz along the c-direction. 
The data for a; = are extracted from Ref.— 



and dielectric constant as function of temperature are 
summarized in Fig. [1]:. It shows that substitution of 
nonmagnetic Ga'^^" for magnetic Mn'^+ leads to reduction 
of T]5^" Hnearly with a slope of -0.6 K per 1% of Ga. This 
is direct evidence of the effective reduction of the strength 
of Mn-Mn interactions, Jmii-Mii and, consequently, we 
suggest that JMn-Tb should be reduced. TpE gets also 
reduced by -1.2 K per 1% of Ga. The most surprising 
effect, however, is that the Ga for Mn substitution affects 
the anomalies at T^^. Indeed, the larger x is the smaller 
and broader becomes the anomaly corresponding to the 
transition to Tb-magnetic ordering (see Fig. [2]). 

Further we studied the magnetic field response of the 
Ga-substituted compounds. Fig. [4] shows isothermal 
magnetization curves as a function of magnetic field for 
all samples. The data on undoped TbMnOa compound 
is in good agreement with previously published results 
For H\\a at 2 K (Fig. [3^), two magnetic transitions are 
observed at i/c 1-7 T and 9 T. The critical fields of 
all mcta-magnetic transitions. He, have been defined by 
the intersection point between two lines: the first be- 
ing a linear fit of M{T) before the anomaly; the second 
passing through the inflection point of the anomaly as 
shown, for example, in (Fig. HJl). The low field meta- 
magnetic transition is a transition to the field forced 
ferromagnetic state associated with the Tb-sublattice. 
Similar meta-magnetic rare earth behavior is reported 
for the isostructural TbAlOs compound»2^ One marks 
out that the magnetization is not saturated up to 14 T 
and the extracted value of the Tb-momcnt, ^ 6.5 jUb, is 
rather low. All these suggest that the Tb moments do 
not become fully ferromagnetically ordered above 2 T. 
Instead, because of their strong Ising anisotropy, they 
presumably arrange in a noncoUincar (Fa, Cy 0) pattern 
(Bertaut notation)^ with rrix ^ 6.5 /ib- Assuming theo- 
retical 9 /iB/Th'^^, that would correspond to 36° canting 
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to a-axis in good agreement with reported data on Tb- 
anisotropyi^Si^ Much less pronounced second transition 
corresponds to the flop of the electric polarization from 
P||c to P||a that is presumably caused by the flop of Mn 
cycloidal plane i^i^S 

Both meta-magnetic transitions are observed in a; = 
4% and x = 10% samples. The field forced Tb- 
ferromagnetic ordering is observed at slightly changed 
critical fields of 1.8 T and 2.4 T for for x = 4% and 
X = 10%, respectively, while the second (polarization 
flop) transition is found at higher field TJ^ = 11 T for 
X = 4% (see Fig. UJ^a) insert). A change in slope is also 
observed for x = 10% around 14 T suggesting that the 
second transition might occur at even higher fields. 

For H\\b at 2 K, also two meta-magnetic transitions are 
observed in TbMnOa in good agreement with published 
datai^ The first one at 1.7 T is a transition to a magnetic 
state where Tb orders with r = 1/3 b*.— The second 
transition is observed at around H\\b ~ 4.9 T. Based on 
the data on TbAlOs^^ we suppose that Tb-spins undergo 
a transition to noncoUinear (C^; F^ 0)-arrangemcnt in this 
state. According to Refi^, this transition is associated 
with the flop of Mn-cycloidal plane from be to ab causing 
the electric polarization flop from P||c to P||ai^i^ Both 
meta-magnetic transitions are observed in x = 4% sam- 
ple. They are shifted to lower critical fields of 1.0 T and 
4.6 T, respectively. Only one clearly separated transition 
with i^c = 4.1 T is observed for x = 10% sample. 

For H\\c^ TbMnOs shows one metamagnctic transition 
as visualized in Fig.lDJc-d) for 2 and 16 K, respectively. It 
corresponds to the transition from cycloidal Mn-magnetic 
structure to simple canted antiferromagnctic structure 
and, simultaneously, to the transition from ferroelectric 
to paraelectric st&te^^ Both substituted compounds 
show similar metamagnctic transitions at lower critical 
fields, Fig.Htc-d). 

In short summary, bulk characterization of the Ga- 
substituted compounds revealed rather similar magneto- 
electric phenomena in all the studied samples. The only 
crucial differences are applied to the Tb-magnctic order- 
ing, where strong suppression of relevant anomalies has 
been observed. To clarify what causes these changes, a 
microscopic view on Mn- and Tb-magnetic orders in the 
Ga-samples has been taken. 



C. Magnetic ordering 

Microscopic insight into the magnetic ordering in 
TbMni_2;Gaa;03 was obtained using ND on powder 
(Fig. [5|) and single crystal samples (Fig. [6]). XRMS has 
been used to contrast Tb-magnetic signal when it con- 
tributed to the same reflections as Mn. Cooling the 
samples below Tj^" we found a set of magnetic satel- 
lites in Brillouin zones ihkl) with extinction conditions 
h + k = even, I = odd, known as A-type reflectionSfSl 
They arise from the ordering of Mn-spins that is incom- 
mensurate (ICM) and characterized by the propagation 
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FIG. 4: Magnetization isotherms of TbMni-ajGai^Oa {x = 0, 
0.04 and 0.1) single crystals with magnetic field along a- (a), 
b- (b) and c-crystallographic (c) directions at T = 2 K. The 
insert in the panel a) shows a zoomed part of H||a magneti- 
zation curve, respectively. The panel d) presents the magne- 
tization isotherm at T = 16 K for H||c. After each hysteresis 
loop the samples have been heated up to 160 K and then zero 
field cooled down to the target temperature. 



vector t'^" = (0 t^" 0). The temperature dependence of 
the propagation vector is shown in Fig. [S^. It decreases 
linearly with decreasing temperature from T^" down to 
TpE for all compositions. Below TpE t stays constant. 
With increasing Ga-content r'^" decreases (Fig. [TJ;) that 
is directly related to change of the Mn-O-Mn bond angle. 
One also observes the reduction of the overall magnetic 
satellite intensities as expected from the reduced total 
magnetization. 

The A-type satellites are typical for all studied com- 
positions in temperature range from Tj^" down to 2 K, 
showing that the main details of the Mn magnetic struc- 
tures for all compounds stay similar. Using the repre- 
sentation analysis f21 the Mn antiferromagnctic ordering 
below T]^" can be described by one single component of 
the irreducible representation r3(0, Ay,0). This gives rise 
to an amplitude modulated AF Mn spin ordering along 
the 6-dircction below T^'^^ 

Below TpE the Mn-magnetic structure of TbMnOa has 
been described by the linear combination of two irre- 
ducible representation Fa © F2 where T2 has also A- 
component along c-axis (0,0,^^)4 The resulting struc- 
ture is a characteristic for ferroelectric phase: cycloid 
?Ti(Mn)=(0, my, mz), where Mn-spins rotate within be- 
plane. The fact that the coupling between Fourier com- 
ponents niz is also of A-type complicates the structure 
analyses as both rUy and rUz contribute to the same A- 
type reflections. However, following a temperature de- 
pendence of a single magnetic reflection often helps to 
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FIG. 5: (Color online) Neutron diffraction patterns of 
TbMni_2,Ga3;03 compounds as measured with E9 diffrac- 
tometer (A =2.816 A) at 2 K. 
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FIG. 6: a) Temperature evolution of the propagation vector 
^Mn q£ TbMui-^Ga^Os with x = 0, 0.04, 0.1 collected with 
the El triple axis spectrometer. Panels b-d) represent the 
temperature evolution of the integrated intensity of G-type (0 
1-T 1), A-type (0 T 1) and A-type (0 2-t 1) incommensurate 
magnetic reflections, respectively. The insert in the panel b) 
zooms the temperature evolution of integrated intensity of 
the G-type reflection around the Tpe. The data have been 
collected on heating. 



see an anomaly at temperature (Tfe in our case) when 
the second magnetic component develops. Indeed, the 
temperature evolution of the integrated intensity of the 
selected (0, 1) and (0, r'^", 1) plotted in Fig.[l:,d 

shows an inflection point at 27.5 K for a; = compound 
in a good agreement with Tfe. Similar considerations 
applied to the Ga-doped samples allow to conclude that 
their Mn-structures have a second magnetic component 
of A-type developing below 25 K and 17 K for x =0.04 



and 0.1 Ga-doped samples, respectively. Both tempera- 
tures are again in perfect agreement with Tfe determined 
from the bulk measurements. The arrows in Fig. ^c) 
show that the relative amount of the spin along the c- 
direction decreases with increasing Ga-content. The re- 
sults of magnetic structure powder refinement assuming 
cycloidal magnetic ordering m(Mn)=(0, niy, rriz) are pre- 
sented in Table [D 

Coupling between Tb- and Mn-magnetic orders man- 
ifests itself below Tj^" where Mn-spin lattice polarizes 
the rare magnetic sublattice^i^ This is reflected in ap- 
pearance of G-, C- and F-type reflections that become 
especially strong below Tfe. Fig. [6)3 shows the tempera- 
ture dependence of selected G-type reflection for all com- 
positions. Unlike the A-type reflections, this one has a 
pronounced concave shape as a characteristic of induced 
order. With increasing Ga-content the intensity of in- 
duced reflections drastically decreases. Since intensity is 
proportional to the squared magnetic moment, ]V[j^(Q)^, 
one can conclude that the effective field acting on the Tb- 
moments from the Mn-spin structure is indeed reduced 
upon Ga-substitution. 

Below T^^ = 7 K, undoped TbMnOs shows very in- 
tense magnetic satellites appearing at positions charac- 
terized by a propagation vector t^^= 0.426(2) b*(~ 3/7) 
and its odd harmonics (Fig. [5]). Quezel et al^ pro- 
posed a sine-wave structure with the Tb-moments lying 
along two symmetrical directions with respect to the b- 
axis (57°) within the a6-plane. According to the results 
from our previous work)^ in one-dimensional case the 
Tb-magnetic structure can be understood as being built 
from |||J.-blocks in the a6-plane with some irregulari- 
ties to adjust Ising-like Tb spins to periodic Mn ordering 
(Fig. [7^). Taking into account data on rare earth behav- 
ior in isostructural 3d-less aluminates RAIO^^^ and 
similar manganites with R = Dy and Ho,— that all 
show the same type of G^Aj, ordering with r = and 
T = 1/2, respectively, one can make this model more 
realistic. Indeed, taking GajAy-coupling to preserve Tb- 
anisotropy and keeping the same periodicity as in Fig. [7^, 
one gets a structure visualized in Fig.[7lD. This structural 
model has been used to fit our powder diffraction data 
and the results of the refinement are presented in Ta- 
ble [11 One has to stress, however, that this solution is 
not unique and the real structure might be even more 
complex. On the other hand, it accounts for both the 
observed harmonic clamping of Tb- and Mn-propagation 
vectors at low temperatures^^ and the Tb-anisotropyi^ 

One of the most striking effects of Ga for Mn substi- 
tution appears below T^^ where magnetic ordering of 
the Tb-sublattice is expected (Fig. [5]). Even the pow- 
der ND pattern of the 4% substituted compound looses 
all the features of long range Tb magnetic ordering: Tb 
superlattice reflections are replaced by low intense and 
broad peaks visible predominantly at low angle positions 
of the strongest Tb-peaks. With increasing Ga-content 
the absolute intensity of Tb peaks drops and their width 
increases (Fig. [T]i). Altogether, the observed behavior 
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FIG. 7: Magnetic unit cell of Tb-sublattice in the a&-plane (in 
the adjacent planes along c-axis the magnetic moments are 
antiparallel): a) The structure proposed by one-dimensional 
ANNNI Ising model from Ref.— ; b) Magnetic structure model 
taking into account Tb-anisotropy 



reflects the destruction of long range Tb-magnetic order- 
ing and its substitution by short range magnetic clusters 
with correlation length less than 25 A (Fig. [TJi). 

The same effect is found by XRMS at the L3 ab- 
sorption edge of Tb. In particular, XRMS data on 
TbMni_j;Ga^03 single crystals showed the presence of 
induced Tb magnetic ordering with the Mn propagation 
vector r'^"=0.262 b* at 2 < T < 28 K for x = 4%. The 
observed magnetic reflections are about 30 times weaker 
that those of pure TbMnOs and exist down to the lowest 
temperatures where no additional reflections correspond- 
ing to Tb long range magnetic ordering (i.e. with t'^^) 
have been observed. For a; = 10%, XRMS found no Tb- 
magnetic reflections in the full temperature range below 
TpE. Single crystal ND, however, reveals some anoma- 
lies in temperature dependence of the Mn propagation 
vector and integrated intensity and width of ICM reflec- 
tions. At 17 K on heating, t suddenly decreases (Fig. [6^) 
and, simultaneously, the width of the all the Mn ICM re- 
flections gets smaller indicating a sudden change of the 
Mn-Mn vs Mn-Tb magnetic correlation. 



IV. 



DISCUSSION 



Based on the results from magnetic and dielectric, 
neutron and x-ray diffraction measurements, a magnetic 
phase diagram of ThMni-ajGa^^Oa compounds with x = 
0, 0.04, 0.1 has been created (Fig. [5]). For small Ga dop- 
ing there is no dramatic change in the magnetic order 
of the system)^ From the dielectric measurements along 
the c-axis we see that Ga doping also does not change 
the ferroelectric state. However, substitution provides a 
platform for studying the interplay between the Mn- and 
Tb-magnetic states. 

According to the phase diagram in Fig. [51 all the mag- 
netic phases persist for x = 4% and 10% compounds. The 
main differences are found in the magnetic state of Tb 
and critical temperatures and fields of the correspond- 
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FIG. 8: Magnetic phase diagram of TbMm-aiGaajOs com- 
pounds 3; = 0, 0.04, 0.1 and magnetic fields applied along 
the H||a (a), H||ti (b) and H||c (c) axes. The squares, circles 
and triangles represent the data obtained by measurements 
of magnetization and dielectric constant for a; = 0, 4 and 
10% samples, respectively. P-M/PFe defines paramagnetic 
and paraelectric phase, while LW-AFM/PE corresponds to 
the phase where the Mn has an antiferromagnetic long wave 
ordering and still paraelectric state. The phases la, lb and 
Ic are identical and correspond to &c-cycloidal Mn-ordering 
accompanied by Ps along the c-direction. Ila is field forced 
ferromagnetic ordering of Tb while lib is an antiferromagnetic 
ordering of Tb with r = i. The Ilia and Illb phases define 
a&-cycloidal Mn-ordering with Ps along a-direction (for x = 0, 
at least). The lie zone is paraelectric phase with simple Mn- 
antiferromagnetic ordering. The dashed lines are expected 
phase boundaries for the Ga-containing compounds. The di- 
electric phase boundary for 1 = has been extracted from 
Ref.i^ 



ing phase transitions. Let us start with Tb-magnetic 
ordering and corresponding effects of Ga for Mn substi- 
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tution. Our results unambiguously show that Jivin-Tb 
is important ingredient for Tb-magnetic order both be- 
low and above T^^. The induced Tb-magnetic ordering 
{t'^^ = r'^" above Tn"*^), should be strongly affected by 
Ga for Mn substitution reducing Jivin-Tb- This finds a 
direct confirmation in our data as both Tfe and the in- 
tensity of the induced Tb-magnetic peaks (i.e. averaged 
induced moment) are significantly reduced upon Ga sub- 
stitution. 

Below T < T^^, when Jxb-Tb competes with JMn-Tb, 
the system minimizes its energy through the matching of 
Tb- and Mn-wave vectors, 3t^^-t^'^ = 1^ The resuU- 
ing Tb-magnetic structure is presented in Fig. [T] Clearly, 
this magnetic regime is supposed to be also strongly af- 
fected by the Ga for Mn substitution. Our data strongly 
support this picture showing that the substitution of 
magnetic by nonmagnetic ions in one of the sublattices 
can change the magnetic ordering in the other sublattice. 
Here, Tb-magnetic ordering becomes of short range in a 
very narrow < a; < 0.04 range of Ga for Mn substitu- 
tion. First of all, such low critical concentration explains 
scatter in the published data on Tb-orderin g^'^^i'^^ since 
small non-stoichiometry in the samples would lead to a 
significant effect on it. Secondly, disappearance of Tb 
long range magnetic ordering proves that Mn exchange 
fields arc involved in Tb magnetic ordering in TbMnOa 
even below T^^. 

Ideally, one can expect by reducing t/Mn— Tb fo get 
an independent Tb-magnetic ordering determined only 
by JTb-Tb- Since TbGaOa does not exist, one can 
again look at Sd-lcss Tb- and DyAlOs both showing 
the same type of GxA.y ordering in ab-plane^^^ Mn- 
containing DyMnos and HoMnOs show the same type of 
magnetic ordering but with doubling of unit cell along 
&-dircctionii^iii Expanding the same considerations on 
TMnOa, one can expect that strong t/Mn— Tb prevents 
Tb of having regular t"^^~1/2 b* structure presumably 
of Ga;Ay-type (in one-dimensional case it would simply 
correspond to ttii)- 

This, however, does not happen since no long range 
Tb-order has been observed at low temperatures. We be- 
lieve the reason for that is Ga for Mn substitution does 
not correspond to homogeneous reduction of Mn-Tb ex- 
change but rather to pronounced effects located at the 
sites which arc occupied by Ga. In this case one has 
to distinguish between Tb ions in the nearest neighbor- 
hood from Ga and not. The former get released from 
Mn-exchange fields, as seen by the significantly reduced 
average induced moment below Tfe, and can form their 
own t"^'^=1/2 b* ordering below 7 K. However, the latter 
follow the Mn periodicity above 7 K and, consequently, 
keep 'Sr^b -r'^" = I'-state below 7 K (Fig. [7b). As a re- 
sult one can still expect a kind of G^Ay (ttii) preserved 
at a local scale while the overall periodicity is completely 
broken. 

This model finds a strong support in the low field data 
on doped TbMui.^Ga^rOs with H||a and H||6 (Fig. 
Indeed, these data show that there are no principal dif- 



ferences to the undoped compound. With H[|a and H||&, 
the low field phase boundary stays almost unaffected (ex- 
cept of H||6 for X = 10% sample where the transition 
becomes very smooth). Having in mind, that this tran- 
sition in TbMnOa corresponds to magnetic field forced 
noncoUinear F^Cj, or C^Fj, arrangement for H||a and 
H||6, correspondingly, one finds it consistent with Tb- 
ordering of G^^A^ (TTii) preserved at local scale in Ga- 
compounds. The high field boundaries for H||a, H||&, 
both getting increased upon Ga substitution, is more 
difficult to interpret. Intuitively, one would expect the 
reduction of both critical fields. However, to answer this 
question one has to know the magnetic and electric states 
of these compounds above these critical fields. This issue 
stays outside the topic of present publication. 

Finally, we summarize the effects of Ga for Mn substi- 
tution on Mn-sublattice. Since it is the Mn-sublattice 
that has been diluted, one can straightforward inter- 
pret the changes directly involving Mn-magnetism. At 
first, one marks out here the reduction of Tj^" and to- 
tal magnetization reflecting the reduced Jmii-Mii- Neu- 
tron diffraction confirms that the magnetic structures ob- 
served below these temperatures are similar while dielec- 
tric measurements show the anomalies pointing to ferro- 
electric character of transitions at Tfe for all samples. 
Among magnetic isotherms one can have a look at the 
H||c data, as Tb-contribution in there is minimal. The 
critical field Hq along c-axis in TbMnOa corresponds to 
the phase transition from ICM Mn-cycloid to a simple 
commensurate antiferromagnetic ordering. Substituting 
nonmagnetic Ga for Mn should result in the reduction of 
He as has been observed in the experiment (see Fig.|3|). 
Generally, one can conclude that all the changes in the 
magnetic properties of Mn-sublattice are mainly of quan- 
titative character and in agreement with effective reduc- 
tion of Jmii-Mh- 



V. CONCLUSIONS 

Ga for Mn substitution in multiferroic TbMnOs has 
been performed in order to investigate experimentally the 
influence of Mn-magnetic ordering on the Tb-magnetic 
sublattice. Single and polycrystalline TbMni_a;Gaj;03 
(x = 0, 0.04, 0.1) compounds have been synthesized 
and characterized by powder and single crystal neu- 
tron diffraction, x-ray resonant magnetic scattering, sin- 
gle crystal magnetization and dielectric measurements. 
The results show that light Ga for Mn substitution does 
not change qualitatively the magnetoelectric properties 
of the Mn-sublattice while it significantly affects the Tb- 
magnetic ordering. The latter looses its long range char- 
acter in a small range of Ga-concentration while pre- 
serves the magnetic structure locally. Thus, our work 
proves that there is a strong influence of Mn-magnetism 
on Tb-magnetic ordering kept down to the lowest temper- 
atures. The whole ordering scheme of Tb in multiferroic 
TbMnOa results from competing Jivin-Tb and JTb-Tb 
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exchange interactions verifying our theoretical model re- 
ported in Refi^ 
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